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Abstract 

We propose a simple analytic representation of the correlation energy ec for the two-dimensional 
electron e; function of the density parameter Vg and the spin polarization C,. This new 

parametrization includes most of the known high- and low- density limits and fits our new fixed- 
node diffusion Monte Carlo simulations, performed for a wide range of electron densities (1 < 
rs ^ 40) and spin-polarization states (0 < C ^ !)• In this way we provide a reliable local-spin- 
density energy functional for two-dimensional systems. The corresponding correlation potential is 
discussed and compared with previous models. 
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I. INTRODUCTION 



The ideal two-dimensional electron gas (2DEG) is a simple model in which N strictly 2D 
electrons are confined in a square of surface S (periodically repeated in space) and interact 
via a 1/r potential within a uniform, rigid neutralizing background. When studying this 
model, one is usually interested in its macroscopic properties, i.e., the thermodynamic limit 
(A^, S ^ oo keeping n = N/ S constant) of its extensive physical quantities per particle. Two 
parameters are enough to define the zero-temperature phase diagram of the 2D electron gas, 
namely the density parameter = Xj^fwaaB (where n is the density and as the Bohr 
radius) and the spin polarization ^ = (n^ — ni)/n, where n^{i) is the density of spin-up 
(down) electrons. 

The model itself is interesting, since it can provide information about electrons confined 
in two dimensions realized in semiconductor heterostructures Q]. Moreover, just like the 
three-dimensional case, the correlation energy of the 2D electron g function of den- 

sity Ts and spin polarization C, provides the local-spin-density (LSD) energy functional for 
density functional calculations of 2D systems. Currently, 2D LSD functionals are based on 
parametrized diffusion Quantum Monte Carlo (DMC) data [2] at C = and Q = 1 [sjO,!^, y|- 
At intermediate spin polarizations, < ( < 1, an exchange-like interpolation is often used 3|. 

We have recently presented new DMC simulations for a wide range of electron densities 
Ts and spin polarizations ( |7i]. The direct DMC calculation of the ( dependence is new and 
provides a reliable basis for building an LSD energy functional for 2D systems. In this work 
we present and discuss an accurate parametrization of these new data as a function of and 
(. This new parametrization accurately reproduces the ( dependence of the DMC data and 
includes most of the known high- and low-density limits. We also compare the corresponding 
correlation potential to previous approximations, finding significant discrepancies at ( ^ 0. 

Hartree atomic units are used throughout this work. 



II. DIFFUSION MONTE CARLO DATA 



which 



Our calculations use standard fixed-node diffusion Monte Carlo (FN-DMC) 
projects the lowest-energy eigenstate $ of the many-body Hamiltonian with the boundary 
condition that $ vanishes at the nodes of a trial function Details of the simulation 
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are similar to those of Ref. P]; further details can be found in Ref. For each of the 
densities corresponding to = 1, 2, 5, 10 we have considered about 20 values of and 10- 
12 polarizations For the densities = 20 and 30 we have used the data of Ref. (sl. We 
have also computed the energy at = 40 for ^ = 1. To estimate the difference A between 
the energy e^irsX) of the finite system and its thermodynamic limit e{rs,() we adopted a 
new strateg y. R ather than a separate size extrapolation for each density based on variational 
energies P.]l^ lll| . we performed a global fit directly based on FN-DMC energies, which 
exploits two physically motivated ingredients: (i) the Fermi-liquid-like size correction 
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A(r„ C, N) = e^irs, C) - e(r., C) = S{1 + AC') [t^(r„ C) - ts(r„ C)] - {v + lC)/N (1) 

{tN and ts are the Fermi energies of and oo particles, respectively, and 6,X,rj,^ are r^- 
dependent parameters); (ii) an analytic expression for e(rs,C), detailed in the next section, 
which involves 12 more free parameters. 

The only uncontrolled source of error, the fixed-node approximation, depends on the 
nodal structure of \E'. We choose a Slater- Jastrow trial function with plane waves (PW) as 
single orbitals. However, within the fixed-node approximation, better results are obtained 



with backflow (BF) correlations in the wave function ll|. Since simulations with the BF 
wave function are considerably more demanding than with PW determinants, we calculated 
BF energies only for ( = 0, N = 58 and ( = 1, N = 57 for each density. For other values of 
and ( the effect of backflow is estimated as a quadratic interpolation in ( and appended 
to PW energies, under the further assumption that the size dependence be the same for BF 
and PW III. 



III. ANALYTIC MODEL FOR THE CORRELATION ENERGY 

In this section we present our parametrization of the correlation energy of the 2D gas 
as a function of and (. We first discuss the ( dependence at a given fixed density, then 
our choice for the dependence is presented, and finally we impose the exact high- and 
low-density limits to our functional form. 
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A. Spin-polarization dependence 



We first noticed that, for Ts ^ 5, the ( dependence of the exchange-correlation energy 
^xc = e — is of our DMC data is accurately described by a biquadratic form, Co(rJ+ci(rs)C^ + 
C2{i^s)C^ (see also Ref. Q]). On the other hand, the known high-density limit 

excirs ^ 0, C) = e,(r„ C) + ao(C) + &o(C)rs Inr, + 0(r,), (2) 

contains non-negligible contributions from higher powers of (: the dominating exchange 
term ex goes like (1 + C)^^^ + (1 ^ O^''^^ and the constant term ao(C) is well fitted by an 
eighth-degree polynomial function of ( [13] . Since we want to interpolate the energy between 
high and low density, we choose a functional form which quenches the contributions to 
beyond fourth order in ( as increases, 

ec(r.C) = {e-^'' - l) (^s, C) + Mrs) + «i(r.)C' + «2(r.)C', (3) 

where 

ef (r„ C) = 6.(r„ C) - (1 + fC' + ifsOe.lrs, 0) 

is the Taylor expansion of ex beyond fourth order in (. Since the first term in the righ-hand- 
side of of Eq. Q contains power 6 and higher of (, it is immediate to identify the function 
tto(^s) as the correlation energy at zero polarization, 

ao{rs) = ec{rs,0). 

Furthermore, 
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2ai{rs) = ^(^c{rs,C) 



gives the spin stiffness, and 



C=o 



B. Density dependence 



We 
Wang 



lave now to fix the dependence of the functions Oj. We generalize the Perdew and 
lJ | form (designed for the three-dimensional gas) to the 2D case as follows 

«,(r,) =Ai + {BiTs + drl + Ar,') In ( 1 + ^ ) . (4) 



This form possesses the small- and large-r^ expansions: 



ai{rs ^0) = A- BiVs Inr, + 0(r,) (5) 

. . . fQ AGA 1 D,F, 1 ... 

a.(r.^oo) = A. + -+^--^j--^-^ + 0(-), (6) 



13l lisj l , provided that the con- 



and it thus has the correct high- and low-density behavior 
straint Ai + Di/Hi = is imposed. 

C. Exact Limits 

Our ec{rs, C) has the correct functional form for small and large Ts] it is now straightfor- 
ward to impose most of the known quantitative constraints. We constrain our e^rs, C) to 
fulfill: (i) the requirement that the exact values Q] of ao(C) ^o(C) at = and 
= 1 in the small-r^ expansion of Eq. are recovered, which implies 

Ao = -0.1925 (7) 
5o = ^(10-37r) (8) 

OTT 

Ao + Ai + A2 + a^PJ^{l) = -0.039075 (9) 

]-3 
127r 



B0 + B1 + B2 = ^^-r^, (10) 



where 



^(C) = (1 + Cf + (1 - Cf' - (2 + IC^ + ^C') (11) 
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(12) 



(ii) the requirement that the total energy e{rs,C) be independent of ( for — »• cxo up to 
order O^rj"^), thus recovering the low-density behavior e —m/rs + njr^J'^ + 0(r7^) 
with positive m and n independent of C^. We thus have 

(13) 
(14) 

(15) 
(16) 









Did 3 
Hi 4""" 


Hx 


C2 


D2G2 3 


H2 


Hi ~ 64""" 
Fi = F2 = 0. 
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We also fixed Ai according to the high-density hmit of the spin susceptibility [13|,|l6|, and 
G2 = since it turned out to be an irrelevant parameter in our fitting procedure. In this way, 
we have built an analytic model which interpolates between the exact high- and low-density 
limits and has 12 free parameters to be fixed by a best fit to our diffusion Monte Carlo data. 
We then perform a global fit (r^, A^), which also includes the infinite size extrapolation of 
Eq. dH), to our data set (122 data for 1 < < 40, < C < 1, and 21 < < 114). In this 
way we fix the values of 36 free parameters, 24 of which disappear from the final analytic 
expression of ec since they only concern the N ^ 00 extrapolation. This fit yields a reduced 
of 3.8. The optimal values of the parameters which yield the model for ec of the infinite 
system are reported in Table B 



IV. LSD CORRELATION POTENTIAL 



The 2D LSD correlation potential /ij!^ for electrons of spin a is given by 

,J(.„0 = 3!!^ . .^(,., _ - (C - =gn.)3^, (17) 

where sgn o" is +1 for spin-f electrons and —1 for spin-| electrons. From our parametrization 
of ec{rs,C) we get: 



decjrsX) ^,, Je-^--(l + /?r,)-l] , ,,2^ ,m 

org 

where J-'{C) and the constant are given by Eqs. (fTTj) and (fT^ . respectively, and 



(18) 



«'.(^s) = = (5. + 2Qrs + 3D,r'g) In 



1 + 



1 



Mrs) 



{B, + C,rl + D,rl)f[{r, 



(19) 

(20) 
(21) 



The derivative w.r.t. C, is simply 



= ^(l-e-^-)r(C) + 2ai(r.)C + 4«2(r.)C^ 



no = i(yr+c-v^ 



c 



2*5 



J_/-3 
16^ ■ 



(22) 
(23) 



It is interesting to compare our correlation potential with the approximations used in previ- 
ous LSD calculations in two dimensions. The most used 2D LSD functional is the one given 



by Tanatar and Ceperley [2], who performed diffusion Monte Carlo simulations at C = 
and C = 1; and gave an analytic fit of the corresponding correlation energies. For the ( 
dependence, many authors used the exchange-like approximation 

ec(r., C) = e.(r„ 0) + - 2 ^f'^^"^^' ~ ^^^^ ^^^^ 

In Fig. n we compare our correlation potential (as a function of and for three different 
values of the spin polarization Q with this widely used Tanatar- Ceperley-plus- exchange-like 
correlation potential. One can see that, while for ( = the two potentials are almost 
indistinguishable, for ( ^ there are significant discrepancies: at ( = 1, the difference 
between the two potentials is ~ 30% at = 1; for lower densities this difference is lower, 
being 15% at r., = 4 and 7% at = 10. At ^ = 1 the discrepancies do not in fact depend on 
the exchange-like choice for the ( dependence: they are exclusively due to the corresponding 
correlation energy of Tanatar and Ceperley, which differs from ours between 35% and 4% 
for rg G [0, 10]. To test the intrinsic quality of the exchange-like interpolation against our 
new (^-interpolation scheme, we plugged into Eq. our new correlation energies at C = 
and C = 1- As shown in Fig. 1 of Ref. j?], the ( dependence of QMC data is rather different, 
especially at lower densities: at the density of the transition to the fully polarized gas, 
Tg ~ 26, the exchange-like interpolation predicts an energy barrier between the C = and 
the C = 1 phases which is more than an order of magnitude higher than the QMC result. 

Our correlation energy which at ( = 1 is, as said, quite different from the Tanatar- 
Ceperley value, should be much closer to the true one since (i) we included the effect 
of backflow on the nodes, (ii) we imposed the exact high-density limit, and (iii) the infinite 
size extrapolation is directly performed on a DMC data set. 



V. SUMMARY AND CONCLUSIONS 



We have presented a new, reliable, LSD functional for 2D systems, based on a new set of 
DMC data for a wide range of electron densities and spin polarizations, and on an analytic 
form which efficiently reproduces these data and includes most of the known high- and 
low-density limits. A comparison of the corresponding correlation potential with previous 
approximations show, for C 7^ 0, differences up to 30% for G [0, 10]. 
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captions 

• Table I - Optimal fit parameters for the correlation energy, as parametrized in Eqs. 
and (jH). Values labelled with * are obtained from exact conditions. The parameter 
Di = —AiHi is not listed (see text). 

• Figure 1 - Correlation potential for spin-up electrons as a function of the density 
parameter and for three different values of the spinpolarization C,. The present 
result is compared with the exchange-like interpolation [3| applied to the Tanatar and 
Ceperley data P]. 
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i = 



i=l i=2 



Ai -0.1925* 0.117331* 0.0234188* 

Bi 0.0863136* -3.394 x 10-^ -0.037093* 

Ci 0.0572384 -7.66765 x 10"^* 0.0163618* 

Ei 1.0022 0.4133 1.424301 

Fi -0.02069 0* 0* 

Gi 0.33997 6.68467 x 10"^ 0* 

Hi 1.747 X 10-2 7.799 x 10"^ 1.163099 

P 1.3386 

TABLE I: 
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